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Abstract 

A new approach to constructing coherent states (CS) and semiclassical states (SS) 
in magnetic-solenoid field is proposed. The main idea is based on the fact that the AB 
solenoid breaks the translational symmetry in the xy-plane, this has a topological effect 
such that there appear two types of trajectories which embrace and do not embrace the 
solenoid. Due to this fact, one has to construct two different kinds of CS/SS, which 
correspond to such trajectories in the semiclassical limit. Following this idea, we con- 
struct CS in two steps, first the instantaneous CS (ICS) and the time dependent CS/SS 
as an evolution of the ICS. The construction is realized for nonrelativistic and relativis- 
tic spinning particles both in (2 + 1)- and (3 + 1)- dimensions and gives a non-trivial 
example of SS/CS for systems with a nonquadratic Hamiltonian. It is stressed that 
CS depending on their parameters (quantum numbers) describe both pure quantum 
and semiclassical states. An analysis is represented that classifies parameters of the 
CS in such respect. Such a classification is used for the semiclassical decompositions of 
various physical quantities. 

Keywords: Aharonov-Bohm effect; solutions of wave equations; uniform magnetic 
field; coherent states; semiclassical decomposition. 

1 Introduction 

Quantum interaction of charged particles with the field of an infinitely long and infinitesi- 
mally thin magnetic solenoid (further Aharonov-Bohm (AB) field) was studied theoretically 
and experimentally already for a long time. In spite of the fact that particle wave func- 
tions vanish on the solenoid line, the particles feel the presence of AB solenoid pQ. This 
phenomenon is called the AB effect and is interpreted as a possibility for locally trivial 
vector potentials to gives rise to observable effects in a nontrivial topology. A number of 
theoretical works and convinced experiments were done to clarify AB effect and to prove its 
existence. By the middle of the 80 's the AB-effect in low energy physics was becoming a 
good instrument for investigating new physical phenomena, principally in condensed matter 
physics, where the AB ring has been the mainstay of mesoscopic physics research since its 
inception, see [2] for a general review. It was discovered that the effect is relevant to a 
number of physical problems, e.g. to anyons in high -T c superconductivity [3], electronic ex- 
citations in graphene with topological defects [H [5] , nanotubes [6] , nonrelativistic scattering 
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in Chern-Simons theory [7J , theory of unparticles [8] , and so on. AB vacuum polarization 
and AB radiation are relevant to cosmic string dynamics, see for example [9| 110] . 

A splitting of Landau levels in a superposition of the AB field and a parallel uniform 
magnetic field gives an example of the AB effect for bound states. In what follows, we 
call such a superposition the magnetic-solenoid field (MSF). Solutions of the Schrodinger 
equation with MSF were first studied in [TT] . Solutions of relativistic wave equations (Klein- 
Gordon and Dirac ones) with MSF were first obtained in [T3] and then used in [T3] to 
study AB effect in cyclotron and synchrotron radiations. On the basis of these solutions, 
Green functions and the problem of the self-adjointness of Dirac Hamiltonians with MSF was 
studied [131 [TSl HS1 [13 EH] • A complete spectral analysis for all the self-adjoint nonrelativistic 
and relativistic Hamiltonians with MSF was performed in [19] . Recently, the interest to MSF 
(and related multivortex examples) has been renewed in connection with planar physics 
problems and quantum Hall effect |20j . It is important to stress that in contrast to the pure 
AB field case, where particles interact with the solenoid for a finite short time, moving in 
MSF the particles interact with solenoid permanently. This opens up more possibilities to 
study such an interaction and correspond a number of real physical situations. For example, 
recent fabrication of a graphene allows one an experimental observation and application of 
effects with relativistic spinning particles under usual laboratory conditions |21) . 

In some cases, it is enough and, moreover, more adequate (and convenient) to use a semi- 
classical description of a physical system. Semiclassical states (SS) of the system provide 
such a description. Usually, such states are identified with different kinds of the so-called co- 
herent states (CS). However, such a formal identification is known for systems with quadratic 
Hamiltonians, in other case both SS and CS construction and their identification is problem- 
atic. In addition to well-known applications of SS/CS in quantum theory [22], there appear 
recently new important applications to quantum computations, see, e.g. [23 . Constructing 
SS/CS for particles in the AB field and in MSF is a nontrivial problem (which was an open 
problem until present), in particular, due to the nonquadratic structure of particle Hamil- 
tonians with such fields. Besides of numerous possible practical applications constructing 
such states gives an important example of SS/CS for nonquadratic Hamiltonians and may, 
for example, answer an important theoretical question: to what extend the AB effect is of a 
pure quantum nature. In a sense constructing SS/CS is a complimentary task to the path 
integral construction, which also is an open problem in the case of the particle in the MSF. 
One can suppose that SS/CS in MSF are in a sense analogous to ones in pure magnetic field. 
In the latter case SS are identified with CS that are well-known, see e.g. [24] [25]. In such CS, 
mean values of particle coordinates move along classical trajectories. The latter trajectories 
are circles whose radii and center position (quantum numbers) label these quantum CS. 
Constructing CS for particles in MSF, we will try to maintain basic properties of already 
known CS for quadratic systems. In particular, such CS have to minimize uncertainty re- 
lations for some physical quantities (e.g. coordinates and momenta) at a fixed time instant 
and means of particle coordinates, calculated with respect to time-dependent CS, have to 
move along the corresponding classical trajectories. In addition, CS have to be labeled by 
quantum numbers that have a direct classical analog, let say by phase-space coordinates. It 
is also desirable for time-dependent CS to maintain their form under the time evolution. 

It should be mentioned that some attempts to construct SS/CS for particles in the 
MSF are presented in the works [26]. However, the states constructed there do not obey 
the principle requirement for SS/CS, the corresponding means do not move along classical 
trajectories. In our recent article [27], we succeeded to construct a principally different kind 
of CS for nonrelativistic spinless particles in the MSF. In the present work we extend this 
construction to the case of relativistic spinless and nonrelativistic and relativistic spinning 
particles both in 3 + 1- and 2 + 1- dimensions (dim.). In addition, developing a semiclassical 
approximation techniques, we constructed SS in MSF on the base of the CS. The progress 
is related to a nontrivial observation that in the problem under consideration there are two 
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kind of SS/CS those which correspond classical trajectories which embrace the solenoid and 
those which do not. It should be stressed that the identification of SS and CS in MSF 
depends essentially on quantum numbers that label these states (on types and positions of 
the corresponding classical trajectories). Particles in constructed SS/CS move along classical 
trajectories, the states maintain their form under the time evolution, and form a complete 
set of functions, which can be useful in semiclassical calculations. In the absence of the 
AB field these states are reduced to the well known in the case of a uniform magnetic field 
Malkin-Man'ko CS [23]. The constructed states give a non-trivial example of SS/CS for 
systems with a nonquadratic Hamiltonian. In addition, they allow one to treat the AB 
effect on the classical language, revealing an influence of AB field on parameters of classical 
trajectories in magnetic field. It should be noted that quantum motion of spin 1/2 Dirac 
fermions is qualitative different in 3 + 1- and in 2 + 1- dim. Since Dirac fermions in 2 + 1- 
dim. (in particular massless ones) describe single-electron dynamics in graphene, we have 
devoted a part of our study to their SS and CS in MSF. 

The article is organized as follows. In sec. [U we start our consideration with classical 
description of particle motion in the MSF. We recall that MSF is a collinear superposition of 
a constant uniform magnetic field of strength B and the AB field (field of an infinitely long 
and infinitesimally thin solenoid with a finite constant internal magnetic flux $). Setting 
the z axis along the AB solenoid, the MSF strength takes the form B = (0, 0, B z ) , 

$ 

B z = B + $<5 (x) S(y) = B H S (r),B = const, $ = const.. (1) 

7rr 

We use the following electromagnetic potentials!]] assigned to MSF ([I]): A = A 3 = 0, 

In Sec. [3J we briefly outline relativistic quantum mechanics of spinning particles in MSF, 
introducing important for our purposes physical quantities. Here we use so-called natural 
self-adjoint extensions of the corresponding Hamiltonians, which correspond to zero-radius 
limit of the regularized case of a finite-radius solenoid. Explicit forms of relativistic and 
nonrelativistic quantum stationary states of spinning particles in MSF are placed in the 
Appendix A. In Sec. IH subsec. 14.11 we build instantaneous CS for nonrelativistic and 
relativistic spinning particles both in (2 + 1)- and (3 + l)-dim., using some universal con- 
structions. In subsec. 14.21 we study the semiclassical approximation. On the base of the CS, 
we construct SS, developing a techniques of semiclassical calculations for means of various 
physical quantities. In Sec. El we construct time-dependent CS for different kinds of parti- 
cles, find and analyze trajectories of means. Some details of these calculations are placed in 
the Appendix B. We summarize and discuss the obtained results in Sec. [6] 



2 Classical motion in MSF 

As was mentioned in the Introduction, our intension is to construct CS in the MSF. Basic 
expected properties of the CS are described in terms of the classical motion in the MSF. 
That is why, we start our exposition with this section, where we present a brief description 
of classical motion of a charge q = ±e with a mass M in the MSF. Trajectories (s) are 
parametrized by the Minkowski interval s and obey the Lorentz equations: 

Mc 2 x u = qF u "x tll (3) 

1 We accept the following notations for four- and three- vectors: a = (a M ,^t = 0, i) = (a , a), 
a = (a 8 , i = l,2,3) = (a 1 = a x , a 2 = a y , a 3 = a z ) , = —a 1 , in particular, for the space-time coordinates: 
x M = (x° = ct, x 1 = x, x 2 = y, x 3 = z), as well as cylindrical coordinates r, ip, in the xy plane, such 
that x = rcosip, y = rsinip, and r 2 = x 2 + y 2 . Besides, dx = dx°d'x., dx = dx 1 dx 2 dx 3 , and Minkowski 
tensor r; M „ = diag (1, —1, —1, —1). 
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where x v = dx^/ds, = d^A v — dyA^. As it follows from ([T]), the only nonzero compo- 
nents of F v ^ are F 21 = —F 12 = B. For trajectories that do not intersect the axis z, we 
obtain from ((3]): 

P = const, P 3 = const, P 1 = exP 2 , P 2 = -exPi, P 2 + P 2 2 = Pi = const, (4) 

where P u = Mcx v — p v — - A" ', is the kinetic momentum, and is the generalized particle 
momentum, P V P V = (Mc) 2 , x — \qB\/Mc 2 , e — sign(qP). Thus, the total particle energy 
£ = cP° is also an integral of motion. One can see that the general solution of (j4|) reads 

. PO P3 . Cp 3 

ct = — — s, z = — — — s + zq = z + zq, x — xq + Kcosw, 

Mc Mc po 

y = yo — ei?sin-0; ^ = xs + ipo = tot + ipo, uj = — — -, (5) 

Po 

where xq, yo, %o, Po, P3, R, and ipo are integration constants. It follows from ([5]) that 
(x - xq) 2 + {y - yo) 2 = R 2 , x = R c cosa, y = R c s'ma , 



2 - x 2 + y 2 = R 2 + R\ + 2RR c cos(ip + eat), R c = \ x 2 + y 2 . (6) 



Projections of particle trajectories on xy-plane are circles of the radii R with central points 
(^o, Uo) placed on the distance R c from the origin. Particle images on the xy-plane are 
rotating with the synchrotron frequency u). For an observer which is placed near the solenoid 
with z > 0, the rotation of the particle with e = 1 is clockwise, and for the particle with 
e = — 1 is anticlockwise. Thus, equations of motion for the charge — q can be obtained from 
equations of motion ([S]) with the charge q by the substitution po by — pq. Along the axis 
z, the particle has a constant velocity dz/dt = —cp^/po. We denote by r max = R + R c the 
maximal possible moving off and by r m i n = \R — R c \ the minimal possible moving off of the 
particle from the z-axis. 

We note that equations of motion for a nonrelativistic particle (P 2 <C (Mc) 2 ) in the 
MSF follow from ([5]) setting po = Mc. Then ct — s and u> = wnr = \qB\/Mc, where wnr is 
the cyclotron frequency. 

The square of the particle rotation energy is E\ = c 2 P 2 L and determines the radius R 
as follows 

P 2 = El (qBy 2 . (7) 
Using (|6]), one can calculate angular momentum projection L z , 

L z = y Pl -xp 2 = ^- (P 2 -P 2 ) + |-, (8) 

which is an (dependent) integral of motion. 

The presence of AB solenoid (the magnetic flux $) breaks the translational symmetry 
in the xy-plane. In classical theory, this fact has only a topological effect; there appear two 
types of trajectories, we label them by an index j = 0, 1 such that j = 1 corresponds to 
(P 2 — P 2 ) > (embraces the solenoid), and j = corresponds to (P 2 — P 2 ) < (does not 
embrace the solenoid), see Fig. [U 

Already in classical theory, it is convenient, to introduce dimensionless complex quantities 
a\ and a 2 (containing the constant h) as follows: 

iPi - eP 2 

ai 



VmM 



ex 




M cx{x - iey) - ■ , . , ,, z , _^ lta 

° 2 = vwrn = Vl/2Rce • (9) 
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i = o 



Figure 1: Two types of trajectories in the MSF 



They define physical quantities R, R c , x, y, and L z as follows: 



R 2 = -a\a\, R\ — 27 1 a* 2 a 2l (x - iey) = V27 1 (a 2 - aj) , 

7 

Pj_ = 2~fh 2 a* 1 ai 1 L z = eh (0X0,2 — atai) + -7 — . 

2nc 



(10) 
(11) 



One can see that aiexp(iwi) and et 2 are complex (dependent) integrals of motion. 

Another important dimensionless integral of motion A (in classical theory A > 0) reads: 



A = 



P0+P3 
Mc 



(12) 



Thus, we can chose the set Xq, yo, zq, A, i?, and"0o as six independent integrals of motion. 
Often, it is convenient to use the light-cone variables x± , 



ct — z, x + = ct + z 



ct 



x + + X- 



(13) 



2 ' 2 
In terms of such variables, the general solution ([5]) takes the form: 

1 + (xR) 2 + A 2 1 + {xRf - A 2 
ct = — 2 x-, z = — x + zo, 

x = xq + Rcosip, y = j/o — tRsinip; tp = ujx_ + tjj , ui = x/X, s = A _1 a;_ , (14) 
where X— plays the role of the time. 



3 Quantum mechanics with MSF 

In quantum theory, it is convenient to represent the magnetic flux $ of the AB solenoid via 
the Dirac's fundamental magnetic flux $0 = ^nch/e as follows: 



($/$o) signB = l a + n 



la 



[($/$o) sigmB] e Z, < n = ($/$ ) signB - l < 1, (15) 



where Zo is an integer and the quantity /1 is called the mantissa of the magnetic flux. In 
fact, /1 determines all the quantum effects in the AB and MSF, see e.g. [15] . We note 
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that the definition (|T5j) differs from the one jl — ($/$o) — [(^/^o)] f° r the mantissa of <£>, 
which was used in some earlier works, and which does not contain the factor signB. The 
quantities /i and jl are related as follows: fi = jl, B > 0; fj, = 1 — jl, B < 0. It turns out that 
the definition (|15j) is very convenient and allows one to write universal expressions for any 
mutual orientations of the uniform magnetic field and the AB flux. 

The quantum behavior of spinning (spin 1/2) relativistic particles in the MSF is described 
by Dirac wave functions $ that obeys the Dirac equation with the electromagnetic potentials 

m, 

ihd t t> = HV, H = c 7 ° (7P + Mc) , (16) 



where 7" = (7°, 7), 7 = (7 fc ) are 7-matrices; P k = p k — ^A k , p k — —ihdk- Below, we 
consider the Dirac equation in (2 + l)-dim., where k = 1,2, and in (3 + l)-dim., where 
k = 1,2,3. 

It is natural that the solutions of the Dirac equation in (2 + l)-dim. and in (3 + l)-dim. 
have much in common. Nevertheless, the algebra of the Dirac 7-matrices for these cases 
is different as well as the spin description, all this implies, e.g., the well known fact that 
quantum mechanics of spinning particle (both nonrelativistic and relativistic) in the presence 
of an uniform magnetic field is essentially different in (2 + l)-dim. and in (3 + l)-dim. That 
is why, we consider the problem under consideration both in (2 + l)-dim. and (3 + l)-dim. 
separately (the former case cannot be extracted from the latter one in a trivial manner). 

It is also known that AB effect in condensed matter physics, in particular, planar physics 
is important in the nonrelativistic case, £±_i a \ ^ Mc 2 . Therefore, we pay a special attention 
to such a limit. In addition, the massless (which is, in a sense, equivalent to the ultrarela- 
tivistic case) Dirac equation in (2 + l)-dim. describes under some conditions the graphene 
physics. In such a case, the Fermi velocity vf ~ c/300 plays the role of the effective velocity 
of light and has to substitute c in all the corresponding expressions. In this connection, we 
consider the ultrarelativistic limit in detail. 

In (2 + l)-dim., a Dirac wave function ^ is a spinor dependent on a; , a; 1 , and x 2 , and 
there are two nonequivalent representations for 7-matrices: 

7 = a 3 , 7 1 = ia 2 , 7 2 = -ia\, C = ±1 , 

where er — (<r l ) are Pauli matrices. Choosing the "polarizations" £ = +1, we describe "spin 
up" particles, and choosing £ = — 1, we describe "spin down" particles. In (2 + l)-dim. these 
are different particles. There exist "spin up" and "spin down" antiparticles. In contrast to 
(3 + l)-dim. case, particles and antiparticles in (2 + l)-dim. have only one spin polarization 
state. 

Stationary states of the Dirac equation with MSF in (2 + l)-dim. have the form: 



= exp 



1>jg(x ± ) , C = ±l, zx = (0,x\z 2 ) , (17) 



where spinors ippj (x±) are subjected to the equations: 



( 



<tP ± + A/or 3 ) il>j»(x±) = Po 4 1 J(x ± ), P ± = (A, A) ; (18) 

(Vo-PV +Mca 3 ) ^(sl) = P o4~ 1) (x ± ) . (19) 

We note that cpo = £ > for particles, and cpo = ~£ < for antiparticle states. 
One can see that 

^ 1) M=o 3 i>^(x ± ). (20) 

That is the reason why we are going to consider only the case C = 1 in what follows. In such 
a case, a self-adjoint Hamiltonian H has the form 

H" = c (<tP± + Mca 3 ) . (21) 



G 



Its domain Dfj depends essentially on the sign d — sign$ = ±1 of the magnetic flux, that 
is why the Hamiltonian has a label §. 

In (3 + l)-dim., a Dirac wave function $ is a bispinor dependent on x , x , x 2 , and x 3 . 
Then unlike (2 + l)-dim. we are not restricted in the choice of evolution parameter of time 
x° , but we can also use the light-cone variable X-. There is also an opportunity to build 
differently spinors adapting them to the nonrclativistic or ultrarelativistic limit. All this is 
described in detail in the Appendix A. 

It should be reminded that all self- adjoint extensions of 2 + 1 - and 3 + 1 - Dirac Hamil- 
tonians in MSF were constructed in jT5J [THl [T5], see, also [T3]. The domains of 3 + 1 - 
Dirac Hamiltonian in MSF are trivial extensions of the corresponding domains mentioned 
in (2 + 1) case, that is why we retain for them the same notation Dfj and we use for 
self-adjoint 3 + 1 - Dirac Hamiltonian the same notation H^. Of course, in this case 
H" 6 = C7° {^2 k=1 2 3 "f k P k + Mcj . In addition, considering a regularized case of a finite- 
radius solenoid, it was demonstrated that zero-radius limit yields two (depending on $) of 
self-adjoint extensions, with domains D^. In contrast to the spinless cas^l both domains 
Dfj involve irregular but still square-integrable radial functions that do not vanish as r — > 0. 
In fact, this means that any wave function is completely determined by its values for r > 0. 
Its value in the point r = can be set arbitrary. In the Appendix A, we represent solutions 
of eqs. (jT5)l and eqs. (TT5|) in (3 + l)-dim. for both values of ■& in the domains Dfj. 

Some important remarks should be made. 

1. In case of spinning particles, some results essentially depend both on the mantissa 
of the magnetic flux ji and on the direction of the flux The latter dependence appears 
due to the spin presence and is specific only for states with irregular radial functions. In 
such states there is a superstrong contact interaction between the magnetic moment of the 
particle and the solenoid flux. Namely this interaction dependent on § can be repulsive or 
attractive. Clear that irregular radial functions appear in the attractive case. 

2. In the spinless case (and even in spinning case) in states with regular radial functions, 
there is a certain translation invariance with respect to a change of the integer number Iq (see 
(|15j) ) by an arbitrary integer k. Such an invariance means that physics depends only on the 
mantissa of the magnetic flux /i. In spinning case, and for fi ^ 0, this invariance turns out to 
be partially broken in states with irregular radial functions (however, a translation invariance 
with respect to the change of Iq by integers k, obeying the condition sign/o =sign(^o + k), 
still holds). 

3. In states with irregular radial functions, we cannot say that particles do not penetrate 
the AB solenoid. However, even in this case, a locally trivial vector potential gives rise to 
observable effects in a nontrivial topology and this is a manifestation of the AB effect for 
such states. 

In (3 + l)-dim., the operators P3, z-component of the total angular momentum operator 

J Z =L Z + £,/2, (22) 

and (aiPi) 2 , where L z = xp y - yp x = -ihd v , otj_ = (a 1 , a 2 , 0), Pj_ = (P 1 , P 2 , 0), 
are self-adjoint on the domain D® H and mutually commuting integrals of motion (all these 
operators commute with the Hamiltonian H ) [TSJ [TB] . In (2 + l)-dim., the total angular 
momentum operator J — —ihd v + ha 3 /2, which is a dimensional reduction of the operator 

J z in (3 + l)-dim., and ^erPj^ are self-adjoint on and mutually commuting integrals 

of motion [HI [16]. One can say that c 2 (crP±^J and J in (2 + l)-dim. and c 2 (a^PjJ 2 and 
J z in (3 + l)-dim. are integrals of motion which play the role of a square of the transverse 

2 The domain of a spinless particle Hamiltonian in MSF involves only regular radial functtions as r — > 0, 
see | 16| . Here we use the terms "regular" and "irregular" as r — > in the following sense. We call a function 
to be regular if it behaves as r c as r — > with c > 0, and irregular if c < 0. 
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kinetic energy {E\ in ([?])) and z-component of angular momentum (L z in ((Sj)) in the case of 
a spinning particle, respectively. Then it is useful to define self-adjoint operators R 2 and R 2 
by analogy with corresponding classical relations (O and ([8]) as follows: 

R 2 = c 2 [aP ± y (qBy 2 , R 2 C - R 2 = -2c (l + ^)h-eJ \qB\~ 1 in 2 + 1 dim., 

i? 2 = c 2 (c* _lP ± ) 2 (gBr 2 , R 2 C -R 2 = -2c Ul +fj,)h- ejJ ^Bf 1 in 3 + 1 dim . (23) 

One can find two types (j = 0,1) of solutions of the Dirac equation which are common 
eigenvectors of operators R 2 and J in (2 + l)-dim. and operators R 2 and J z in (3 + l)-dim., 
see (|97|) , (|130p . and (|136[) in the Appendix A, respectively. Such solutions have two quantum 
numbers n\ and n,2 in common then we may be using to them the general notation $fn' M2 . 
Note that eigenvalues of the operators J and J 2 are the same, J z = J = eh(lo — I + 1/2), 
where I is an integer. Then, using an appropriate inner product on ccy-plane, see (I108|) in 
(2 + l)-dim. and (|127[) in (3 + l)-dim., we obtain the mean of the operator R 2 — R 2 , 

(*,*) = (*, in2 + l dim., = in 3 + 1 dim . (24) 

In the semiclassical limit the sign of the mean allows one to interpret the corresponding 
states as particle trajectories that embrace and do not embrace the solenoid. Namely, an 
orbit embraces the solenoid for I > (type j = 1) and do not for / ^ — 1 (type j — 0). This 
classification corresponds to the classical one introduced in the previous section, see eq. ([8]l 
and Fig. [TJ Trajectories with / = 0,-1 are situated most close to the solenoid. 

For ji = there is no any impact of AB solenoid on the energy spectrum and the energy 
spectrum is given by the Landau formula. For /i ^ 0, energies of states with j = 1 differ 
from the Landau levels, whereas energies with j — coincide with the Landau levels, see 
(|107p in the Appendix A. 



4 Instantaneous CS on xy— plane 
4.1 Quantum states 

In spite of the differences between stationary states of spinning particle in (2 + l)-dim. and 
in (3 + l)-dim., see (1M1) . (jl33[) . and (|135j) . in both dimensions one can build CS on icy-plane 
in a similar manner, using some universal constructions. Let us introduce operators ax, ci2, 
and a\,a\ that correspond to classical quantities ai, a2, and a*, a\ , 

iP x - eP 2 _ Mc x(x- iey) - iP x - eP 2 
^2hMcx ^ ~ V2hMc>c 
iPi + eP2 .t = Mcx(x + iey) + iP 1 -eP 2 
\f2hMcK ° 2 V2hMcx ' ° 

It should be noted that operators Pi and Pi are symmetric but not self-adjoint on the 
domain D^. That is why, one cannot consider a\ and a\ as adjoint to di and a-i respectively. 
Nevertheless, the operators (|25[) play an important role in the further constructions. 

Using properties of Laguerre functions, one can find the action of these operators on the 
functions fllOip : 

ai^?,„ 2 ,a(¥>, P) = y/niQnl-iwfa P) ' "'t^nln^A'P' P) = Vni + l^+i,^,,^, p) , 

(26) 



ai = 
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where possible values of n\ and n 2 depend on m,l,a, and j according to (|101[) and the 



functions $ 



(j) 

ni+si, ri2+S2,<7 



are defined in p04p . 



Formal commutators between the operators a\,di, and a\ . a 2 have the form: 



01, o{ 



1 + /, 



a 2 , a 2 



1-/, [ai,a 2 ]=/, 



ai, a 2 



0. 



with a singular function / = $ (7r_Br) _1 <5(r) = 2(/o + fi)S(p). All the solutions of the Dirac 
equation in (2 + 1)- and (3 + l)-dim. introduced in the Appendix A are expressed via the 
functions (|101|) and describe states of spinning particles out of the solenoid, r > 0. In such 
a case, the function / gives zero contributions and can be neglected. In turn, this means 
that on the domains Dfj operators a\,a\, and 0,1,0,2 behave as creation and annihilation 
operators. Then quantities x — iey and L z can be expressed in terms of the operators a\, Si 
and a\, 0,2 as follows: 

(a; - iey) = (a 2 - afj , ~L 2 + l Q + f i = (n x - N 2 ) , N s = 0%, s = 1,2, (27) 

which is the same in 2 + 1 and 3 + 1 dim. The operator c 2 (^crP± S j in 2 + 1 dim. case and 
c 2 (a^PjJ 2 in 3 + 1 dim. case have the following representations: 



c 2 (o-P_l) = 2/ic \qB\ TVi + (1 - a 3 e) /2 
c 2 (a ± P ± ) 2 = 2hc \qB\ \Ni + (1 - S z e) /2 



in 2 + 1 dim . 
in 3 + 1 dim . 



(28) 



where relations (f26| and expression (|107|) are used. Then the operators R 2 and i? 2 have the 
following form (see the definitions (f2"3|) ): 

fi 2 = 7- 1 ^2iVi + 1 - o- 3 e) in 2 + 1 dim., ^ 2 = 7" 1 ^2iVi + 1 - £ z e) in 3 + 1 dim., 
R 2 C = 7" 1 (27V 2 + 1^) in 2 + 1 and in 3 + 1 dim . (29) 

As was already mentioned in the Introduction, our aim is to construct CS. One can for- 
mulate a definition of CS for systems with quadratic Hamiltonians, see [55]. Unfortunately, 
no general definition of CS for arbitrary quantum system exists. In our case, with a non- 
quadratic Hamiltonian, defining CS, we would like to maintain basic properties of already 
known CS for quadratic systems. First of all, these states have to minimize uncertainty rela- 
tions for some physical quantities (e.g. coordinates and momenta) at any fixed time instant. 
Second, means of these quantities, calculated with respect to time-dependent CS, have to 
move along classical trajectories. It is also desirable for time-dependent CS to maintain 
their form under the time evolution. Here, it is supposed that time-dependent CS are solu- 
tions of the corresponding wave equation, Dirac or Pauli (and Klein-Gordon or Schrodinger 
equation in the case of spinless particle). Thus, the problem of constructing the CS states 
is in main reduced to a suitable choice of the form of the CS at a fixed time instant. We 
call such CS instantaneous CS (ICS) in what follows. In the case of quadratic systems, e.g. 
a non-relativistic particle in the magnetic field, ICS are eigenvectors of the corresponding 
annihilation operators, let say the operators a±, a 2 from (|25[) without the magnetic flux. 
In the case of MSF where particle Hamiltonians are nonquadratic, these operators are not 
exactly annihilation operators for both types of the functions $nZn 2 ,<7- Nevertheless, as 
is demonstrated below, one can construct some kind of ICS that have the above described 
properties. These states are very close to eigenvectors of the introduced operators ai,a 2 
from (1251) . At the same time, these states maintain their form under the time evolution. 
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The presence of the AB flux breaks the translational symmetry in the xy-plane. That 
is why in the problem under consideration, there appear two types of CS, those which 
correspond classical trajectories which embrace the solenoid and those which do not. Taking 
into account classification of quantum states according to types j = 1 and j = 0, which 
depends on the sign of the mean value ([24")) , we see that each of these CS must be constructed 
using stationary states of the same type. 

It is convenient to pass from the functions <&nl,n 2 ,o IjlOlj) to new functions 3>zP )Z2)CT as 
follows: 



Z?Z?*$, na A<P, P) 

v /r(i + ni)r(i + n 2 )' 



(30) 



Here z\ and z 2 are complex parameters, possible values of n\ and n 2 depend on m, I, a, and 

j according to eq. (jlOljl . and we set Af = 1. The functions 3>aj,* 2 ,<7 can be expressed via 
special functions Y a , 

*^U,«r(Vi (°) = ex P i ze Po — i + (1 — e^) /2] ^} F_ Q? (zi, z 2 , p) , aj = I- (1 - ecr) /2 + /i, 
*S r »,«r(v. P) = exp [h - I + (1 - «x) /2] p + 7T [I - (1 - e) (1 + a) /4]} F ttf (z 2 , z 1; p) , 

~m „m+c* r / \ 

6\ ^ 2 ^m+a,m\rj 



z 2 ; p) = > 

^ oV /r(i + m)r(i + m + a ) 

By the help of the well-known sum, 

■2* ^a+m,m(^') 



(31) 



E 



r(, v/r(i + m)r(l + a + m) 



= z 2 cxp (z — x/2) J a (2y/xz ), 



where J a are the Bessel functions of the first kind, one can obtain the following representation 
for Y a : 

Y a (zi,Z2jp) = exp (ziZ2 - ~J (^/zz/zij J a {2 v /ziz 2 p). (32) 

Using the functions <&ii, Z2 ,o-, one can construct ICS on xy— plane. In (2 + l)-dim. and 
in (3 + l)-dim., ICS are constructed by the help of spinors described in the Appendix A by 
substituting the functions $^1,712,0- for the function Z2i(T . 

Thus, using eqs. (|55]l. (|100l) . and ([55]). we obtain ICS for massive, ( = +1, spinning 
particles on xy-plane and in (2 + l)-dim.: 

P) - {* 3 [±Ho (#1) -o-Px] + Mc}ug ) ifti±1 fop), 



TVi ) = fig— / e - n o- 2 dr, fig = M 2 c 2 + 2h \qB\ /c JVi + (1 - ere) /2 
7 V 7 *" Jo L 



For such states, we have 



= 2Mc * 



2 / ZJ 



■ 2i,Z 2 „±ll 



n JVi + Mc 



0') A 



(33) 



(34) 



where the inner product is defined by (|105[) . 

According to (|112j) . (| 1 13[) . and (| 1 14|) . nonrelativistic ICS for 2 + 1 spin up particles have 
the form: 



(35) 
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whereas for the spin down particles they read: 



( mW down vi>(j') dowl 



2 / D 



■ zi,z 2 „Tl' z'^z'^l) j_ ' 



(36) 



According to (|138[) . ICS on xy— plane for 3 + 1 nonrelativistic particles with a given spin 
polarization s = ±1, have the form: 









, T ,(i)down / \ h 

W P)/ 



XR 



y zi,Z2,, 



±s) zj,z 2 ,±s 



(37) 



where the inner product of four-component spinors ^ and ^' on xy-planc is defined in (|127|) . 
According to (|115[) - (|117[) . ICS for 2 + 1 massless £ = +1 fcrmions are 



V+«e*kU,-i(¥>> P) 



*£!, 2 , + i(^p) 



,(0, + l) 



^"Lefo P) 



SLZa.+lWPj \ (o.+l) / \ 
= (0) , J + U 0,zJ 2 WP), 



V /r ( 1 +"2) 



2 v /r(l + n 2 J 



m=0 



The inner products of such states have the form 

y ±,zi,z 2 > ±,z[,z' 2 J D \ Zi,Z2„v' z[,z 2 ,aj ± 



(38) 



(39) 



a=±l 



In the same manner, by the help of (|118[) , one can obtain ICS for 2 + 1 massless ( = — 1 
fcrmions. 

According to eq. (|133p . for 3 + 1 relativistic spinning particles, ICS on xy-plane have the 
form 

UWJ^fo P ) = *E 2 >, P) ( _X ) ■ (4°) 
The inner product of such states on xy-plane reads: 

One can see that in all the cases, the inner product of ICS on xy-plane is expressed via 
the matrix elements: 



= {sfWn V^R) , Pa = P - (1 - 0a) ; 

oc 

Q a {u, v) = Q~(u, v) + (v/u) a I a {2uv), Q-(u, v) = ^ («/< + ' 4 +! (2w), (42) 



z=i 
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where I a are the modified Bessel functions of the first kind. We note that in contrast to 
the spinless case [27], subindex a in the functions Q a (u, v) can take also negative values 
-Ka<0. 

It should be noted the importance of the obtained result. It turns out that all the means 
and matrix elements with respect to the ICS are expressed only via two functions Q a {u 1 v) 
and Q~ (u, v). That is why the further study of such physical quantities is reduced to the 
analysis of these functions. 

It follows from (gU) that: 



(43) 



and 



a 2 ^\ Z2 J^ p) = z 2 *£\ Z2 J<p, p) + (-l) j #& I(X (V, P) 
Eqs. allow one to calculate the matrix elements 



i = -(l+0€)/2 

r=(l-tfe)/2 



aL<f> U ' ] , 



(44) 



(45) 



Results of such calculations are, for example: 

(46) 

Using an appropriate inner product, see above, we define means of an operator F with 
respect of the ICS on xy-plane, (F)^y Then, we consider important cases when a matrix 

operator F is either the identity matrix / multiplied by a differential operator f,F = fl, or 

F = c 2 (crP±j in 2 + 1 dim., and F — c 2 (a±Pj_) 2 in 3 + 1 dim. Here, we can express (F)^ 

via the means \ $*i,z2,<t, i^-**?^,^ , where either F a = /, or F a = c 2 P\ — ehc\qB\a. 

Thus, we obtain, for example, for 3 + 1 relativistic and 2 + 1 nonrelativistic spin up particles 
the following expression (with the corresponding interpretations of the number a) 



I gSj) *0') 

1 ^Zl,Z2,CT) *Zl,Z 2 ,ff 



(47) 



For 2 + 1 massive £ — +1 relativistic particles we obtain 



(*& a ,a, [n (^<) + Mc] *g> 



(48) 



whereas for 2 + 1 massless ( = +1 fermions, we have: 



6&zi, Z 2 ,<7 



U) 



(J) 



(49) 
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Note that means (F) for 3 + 1 nonrelativistic spinning particles and antiparticles at given 
s are expressed via means (|47|) for 2 + 1 nonrelativistic case according to ([37]) . 

Then, using (j43|) and notation (|47|) . we obtain the means of operators Nk, k = 1,2, for 
example, 

(ivfe)^ =%a^ln^ (50) 



for 3+1 particle and for nonrelativistic 2 + 1 spin up particles; 

(iVfc) w = z k d z . In J] fcCO 



cr=±l 



for 2 + 1 massless £ = +1 fermions; 

2fe 9 4 ln{[n (49 zl )+Mc]^«} 



for 2 + 1 relativistic massive £ = +1 particles. 
Using (HU), we nn d that 



(x - iey) y) = y/2j 1 (02) m - (01) 



C?) 



(51) 



(52) 



(53) 



where, for example, 



(°i)(o) = 2i A l-/v(M>N), (a 2 ) m = z 2l 



(0) 



Q-(u v) 

(ai) {1) =zi, (a 2 ) (1) = z 2 A^(\z 2 \, \zi\), A a {u,v) = a 



(54) 



for 3 + 1 particle and for nonrelativistic 2 + 1 spin up particle (a = +1) and antiparticle 
(a = —1), and 



7-^ _ E a =±iQi-^(\ z ^\^\) --^ 



(55) 



for 2 + 1 massless £ = +1 fermions. 

Note that one can get the means {Nk)^ and (ttfc)(j) f° r the case of spinless particle from 
expression ([SO]). (|54]l . and ((42]) at = see [27]. 



4.2 Semiclassical approximation 

Representations ([5"U|) - (T52~j) allow us to connect means of i? 2 and -R 2 with the parameters z\ 
and z 2 . It follows from (f2l?)) that 



2(JVi) 



0') 



1 - ere 



(56) 



Note that these relations are valid in the case of spinless particle at a — 0. We expect that 
in the semiclassical limit iNk)u\ ~ l z fc| 2 - At the same time length scales defined by means of 

(R 2 )rj) and (Rc)(j) have to be sufficiently large, which implies \zk\ 3> 1 in the semiclassical 
limit. 

We note that in the pure quantum case, as a characteristic quantum scale of the rotational 
motion we can take the quantity 



£^ uant = 2 \qB\ he = 2MV \B\ /B , B q = M z c A / \q\ ft, 
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where Bq — m^c 3 /eh ~ 4, 4- 10 13 G is the critical magnetic held above which the nonlinearity 
of QED becomes actual. The corresponding length scale is 



Equant = V27^ = y/2B /\B\X c , Ac = h/Mc. (57) 

For the angular momentum projection J z a characteristic quantum scale is obviously h. For 
a given energy, i.e., for a given (i? 2 )^, the quantity {Jz)(j) is proportional to {Rc)(j) ^ ue 
to (|23|) . and, therefore, can be characterized by the corresponding length scale i? qua nt- Note 
that the i? qU ant is much larger than the Compton length Ac if the magnetic field B is weak, 
B / \B\ ^> 1. Thus, the conditions \zk\ 2 ^> 1, correspond to ones 



quant 



At the same time, in the quantum case, the dimensionless quantities \zk\ are of the order 1. 
We see that the semiclassical decompositions are adequate namely in case of strong enough 
magnetic fields (e.g., pulsar magnetic fields B for which Bq/ \B\ ~ 10 -2 and i? quan t <C Ac)- 
We expect that the sign of the difference 



(58) 



is related to the trajectory type in the classical limit. One can see that such a limit implies 
the following conditions: 



Hi) 



> R, 



quant 



M|»i. 



In particular, for states with j = 0, we have \z\\ <C \z2\, and for states with j = 1, we have 
3> \ z 2\- We note that in both cases the corresponding functions Q a (u,v) are calculated 
at M > lul > 1. 



There exist all the derivatives d v 



and the series YlbLi ®v (v/u) a+l I a +i(2uv) 



(v/u) a+l I a+ i(2uv) , the series Q a (u, v) (|4"121 converges 



converges uniformly on the half-line < Re v < 



oo. Thus, one can write a differential equation for Q a (u, v), 

d v Q~(u, v) = 2v [(v/u) a I a (2uv) + Q~(u, v)] 
A solution of this equation, which corresponds to (J42J , reads 



Q~{u, v) = 2e v / e" 5 {v/uf I a (2uv)vdv . 
Jo 

Using eq. (6.631.4) [28], we represent this solution as follows: 

Q~(u, v) = e u2+v2 Q-(u, v), Q-(u, v) = [1 - T(u, »)] 



T(u, v) = 2e 



e v (v/u) a I a (2uv)vdv. 



(59) 



Then 



) a (u, v) 



t (u,v), Q a (u, v) = l — T(u,v) + e 



—u —v 



(v/u) a I a (2uv) , (60) 



and we can calculate, for example, the means (|50[) 

2 



= Qi-n. > = (a/44, v^ii) 



(61) 
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These means can be represented explicitly in an real form, taking into account that 

4,i l-gxj d u Qi-^(u, v) + 4,2 \z2\dyQi^Ju, v) 



z k d z , k lnftW 



2Qi-^(u, v) 
4,i \zi\ dyQ^tu, v) + 4,2 \zz\ duQ^iu, v) 



u=\z 1 \, v = \z 2 \ 



2Q MtT (u, v) 



u=\z 2 \, v=\zi\ 

0. Thus, the contribution due to the 



We stress that means (|5T|) allow the limit 
AB field can be easily isolated. 

Using power decomposition of the function v a+1 I a (2uv)e~~ 2uv near the point v = v for 
an estimation of the integral T(u, v) in (|59p and asymptotics of the function I a (2uv), one 

for It; I > Id 3> 1. Thus, we obtain the semiclassical 



can see that Izkf 3> z^ d z > Inila^ 



expansions: 



\Z2\ 



\Zk\ » 1> 



(62) 



which connect means of i? 2 and i? 2 with the parameters zi and 22- Thus if ||zi| — \z%\\ 3> 1, 
eqs. (1621) imply the following relations 



0; \z 2 \, j = l. 



(63) 



It should be noted that relations (|63|) have nothing to do with conditions of the ap- 
plicability of the semiclassical expansions (|62j) . Obtaining the latter expansions we have 
supposed that \z\\ < |z 2 | for states with j = 0, and \z\\ > \z%\ for states with j = 1. There- 
fore, relations (|62p between means of R and R 2 and parameters Z\ and 22 take place even 
if a definite relation between sign ( \z\ \ — \z 2 \) and j is absent. 

Retaining only leading terms in decompositions (|62|) . we reproduce the corresponding 
classical relations (fTOj) with \z\\ — \a\\ and \z 2 \ — |a 2 |- In other words, one can say that the 
classical relations (|10[) correspond to the leading approximation for sufficiently large radii. 
Thus, the leading approximation in the semiclassical expansions corresponds to the classical 
limit. Next-to-leading terms define physical quantities in the semiclassical approximation. 
These terms depend on the space dimension and particle spin. 

Let us consider semiclassical approximation retaining next-to-leading and next-next-to- 
leading terms. If \v\ 3> \u\ 1, one can approximate the integral T(u, v) in (|59j) by a power 
series in u/v as follows: 



T(u, v) = (v/uf I a (2uv)e~ u ~ v (1 + u/v + 
Then, using asymptotic of I a (2uv), we obtain from (|60|) : 



•)• 



Q a (u, v) = 1 



(u/v) 



1-a 



which implies 

d v Q a (u, v) 



d u Q a {u, v) (u/v) 



1/2- 



-(«-• it) 2 



-{v- 



\v\ > Id > 1. 



(64) 



(65) 



Qa(u, V) Qa(u, v) \fK 

Thus, for semiclassical states corresponding to orbits situated far enough from the solenoid, 
i.e., for \\zi\ — \z 2 \\ ^> 1, the terms Zk d z > lnT^r are small as exp (— |zi| 2 — \z 2 \ 2 )■ 

Then the semiclassical expansions (|62p in the next-to-leading approximation reads: 
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(R 2 ) U) -(l-ae)/2, 



\Z2\ 



= l(^)(i) 



1/2, 



zi 



Z 2 >1. 



15 



In the most interesting case when a semiclassical orbit is situated near the solenoid, such 
that the condition ||u| — \u\\ <C 1 holds, the influence of AB solenoid (due to // ^ 0) on the 
orbits are not small. In such a case 

~ 1 v — u a — 1/2 

Qa{U, V) « - + —j= , 

d v Q a (u,v) ^ duQgju, v) ^ 2 / 2 ^~ M I Q ~ 1 / 2 ^ (66) 
Q Q (u, u) ~ Q a (u, v) ~ V >/™ / ' 

such that, for example, means (|6ip are 



(jvjo W) « i«*r + (-i) 



(-iy +2 



2tt 



|zi| - |z 2 || « 1. (67) 



Thus, for \ \zi\ — l^aH -C 1 all /i-dependent contributions to the means (Nk)^ are of the order 
one, which is natural for a pure quantum case. Next-to-leading contributions to the means 
(of order \zk\) that does not depend on fj, are much bigger. These semiclassical contributions 
appear since each of j-type ICS includes only a half of eigenfunctions of the operator J z . It 
follows from (|67[) that in the leading approximation 



(iVi) m - (JV 3 ) m « (-l) i+1 |Zll t- |Z21 , IN - N2II « I- (68) 



'0') _ 



in 



At the same time, relations (|62l) yield in the semiclassical approximation: 



kii 2 «i(fl%)+(-i) , 7^(^) 0) , 

M 2 « - (-i) J IN - Nl « l. 

Then, using ([68]) . we obtain 

IT 



- - {V( R2 h) + \/( R2 ch)) > 11*1 - 1*11 « L 

which implies 

- £ « * (-iy 4 (x/wh, + JrmQ Ku\ nt . (69) 



~ 2^ ~ v ; V V (i) V W / quant 
Then the quantity d(j) (|58f is: 

d {j) *(-l) j+1 ^, ||z 1 |-|2 2 ||«l. (70) 

Thus, for ||zi| — 1 1 <C 1, and in the semiclassical approximation, the mean minimal 
possible moving off |<i(j)| of the particle from solenoid line is of order i? qU ant, in particular, 
d(j) < for states with j = 0, and drj) > for states with j ' = 1, independently on the sign 
of the difference N — \z2\. 

Eqs. (|43)) and ([60]) allow us to calculate variances of the operators Nk, 



Var 3 
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In the semiclassical approximation, we have 



Var, (N k ) « \z k \ 2 , \\ Zl \ - \z 2 \\ » 1; Var, {N k ) « (1 - 1/tt) \z k \ 2 , \\ Zl \ - \z 2 \\ « 1. (71) 

Thus, standard deviations of R 2 and R 2 in the semiclassical ICS are of the same order 
for any value ||zi| — \z 2 \\, namely: 



Si (i? 2 ) - ^Var, (R 2 ) ~ iwft ^ (#c) = y Var, (i??) ~ i? quant 



C/O- 



In this case, the typical spread of the radii i? and i? c are given by the standard deviations 



5 j {R) = S j (If) (R 2 



'(j) 



-1/2 



ft 



quant 7 



^ (i? c ) = ^ (i? c 2 ) {R 2 C ) U) 



-1/2 



I?. 



quant ■ 



(72) 



For — |z 2 || <C 1, the difference (R 2 )^ — (R 2 )^ is of the order of the standard 
deviation of R 2 —R 2 , which is Sj (R 2 ) +8j (R 2 ) ■ Therefore, the mean angular momentum (IMf 
is of the order of J z standard deviation, and for \ \zi \ — \z 2 \\ -C 1 , the quantum scale of angular 
momentum is much greater than h. In this case \ d^) | is of the order 5j (R)+Sj (R c ) ~ -Rquant- 

We note that ICS a la Malkin-Man'ko [24] (the case $ = 0) can be associated with 
the superposition of the j = 1 and J = states for /z = 0, which includes already all the 
eigenfunctions of J z . The inner product on the xy-plane between such states is the sum 
K ( v ] +ni 1) at fi a = 0. Using eqs. (8.511.1) from we hnd 



Therefore, in such ICS (N k ) = \z k \ . Similar mutual compensations take place in expressions 
for 7?i ^and for fi ^ in the semiclassical limit, \z k \ 2 ^> 1. E.g., for ||zi| — |z2|| -C 1, 
we obtain: 



n (o) +n (i) = 1 + 



-Si- 



One can see that in states that are superpositions between different j, leading corrections 
to means (N k ) = \z k \ 2 (|67| disappear in the semiclassical approximation. 

If \zi | and | z 2 1 differ essentially, i.e. , 1 1 z\ \ — \z 2 \\ 1 , one may believe that next-to- leading 
terms in 1Z„ \ given by (|64|) . remain uncompensated. That is not true. To see this, one has 
to take into account that next-to-leading terms in TZ„ + TZ„ are due to contributions from 
(ffi4"|) and from leading terms in Q a (u, v) for \u\ 3> \v\, 



) a (u, v) 



1 



2J-KUV 



(v/u) a e 



-(v-u)' 



, M > |u| > 1. 



(73) 



We recall that for ICS, the domain \u\ > \v\ is not classical one even if \z k \ 3> 1. 
Using (l53l) . (J54J) , and the representation 



A Q (u, v) = 1 — zu a (u, v), w a (u,v) = (v/u) a I a (2uv)e 



Q a (u, v), 



we find: 



(x - iey) U) = [(a 2 ) w - (ai) (j) 

(02) (0) ~ («i)(o) = z 2- z*Ai-p(\zi\, \z 2 \) = z 2 - z{+ zlzui-^(\zi\, \z 2 \ 
(a 2 ) (1) - (ai) (1) = z 2 A fi (\z 2 \, \zi\) - z* = z 2 - z* - z 2 w^{\z 2 \, \zi\). 



(74) 
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By the help of eqs. (|54|) and (l67l) we calculate the variance of (x + y) in j-type states: 

2 



Var,- (x + y)= [\x- iey\ 2 ) - (x- iey) 
v 7 u) 



2 7 - 



+ (iV 2 ) f 7l + 1 - (a x ) 



a) 



(a 2 ) 



0) 



(75) 



Let us consider the semiclassical limit |z/J 2 ^> 1 for ICS with j = and < 1^1 and for 
ICS with j = 1 and |zi| > |z2| , in both such cases \v\ > |it| 3> 1. In this case, using the above 



M 0) - (ai) 



(i) 



Z2-Z1 



results, one can verify that corrections to the classical expression 

are small. In particular, using eqs. (|64|) and (1661) . and asymptotics of I a (2uv), in the next- 
to-leading approximation we obtain the following result 



d a (u,v) 
d a (u,v) 




u\ < 1. 



Thus, eqs. (|74j) match with once ([6T| in the classical limit, and due to (J61 
in the semiclassical approximation, the variances (|75[) are relatively small, 

Varj (x + y)Ks 2 7 ~\ \z 2 \ > M for j = 0, |zi| > \z 2 \ for j = 1; 



(76) 

S3, and (l66j) . 



Var (x + y) 



4|*i| 



Vari (x + y) 



. IN 



Ml < 1. 



(77) 



' %/^7 ' 

However, near the AB solenoid, where ||zi| — \z%\\ -C 1, the variances increase significantly. 
Thus, it turns out that quantum length scale for *J (i? 2 )^ and \J (R 2 )^ is essentially 
different from quantum length scale for (x)y\ and (y)^y 

Note that in this case, the principal part of next-to-leading contributions to 0*0 (j) and 
(y)(j)> gi ven by ([76]) . do not depend on /i, which is quite similar to the behavior of (-Wfc)(j) 
given by (|67|) . However, here a continuous limit to the case fi = does not exist. This can 
be checked considering means (x) and (y) in the superposition of the j = 1 and j = states. 
For example, the mean (x) in the latter superposition includes both means (x) ,^ and some 
interference terms given by (|45j) . The latter terms are absent only for fj, = 0. That is why 
the means (x) and (y) in the Malkin-Man'ko CS cannot be obtained in the limit /x — >• 0. 
Thus, namely means {x)rj\, (u)(j)i an d (x), (y) are especially sensitive to the topological 
effect of breaking of the translational symmetry in the xy-plane due to the presence of AB 
solenoid. 

In the pure quantum case, the mean values depend significantly on the particle spin 
and on the mantissa /i and are quite different from the corresponding classical values. E.g., 
for small \u\ and \v\ the functions Q a (u, v) with positive and negative a behave essentially 
different. Using the representation (8.445) [28] for the function I a (2uv), we obtain: 



t (u, v) 



^ 2a /r(l + a), |w| < 1 

. . |u«|-K) 

> a {u, v) — > 0, a > 0, 

|u«|-HI 



I Qa(u, v) ' QO, a < 0. 

Note that in the case of spinless particle a > (see [27 ), while a can take also negative 
values —1 < a < in the case of spinning particle. Thus, e.g., for [ziZz] "C 1, ( a 2)(j) ~ ( a i)(j) 



differs essentially from z 2 

( a 2) (0 ) 



i . We have A a (u, v) = v 2 (a + 1) . Then, 
(ai) (0 ) ~ 22, (a 2 ) (1) - (ai) ( i) « 



(78) 
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In this case, using (|43|) and (|42|) . we obtain Varj (x + y) sa 27 _1 . However, it is big in 
comparison with small \zk\ 2 ■ 

If |u| 3> \v\, we deal with the quantum case even for big \zk\ ■ Here A. a (u, v) = v/u — > 
which gives a justification for relations (|78p . In addition, in the quantum case, we have 
(Nk)(j\ ~ 1, and, at the same time, contributions to C^fc)(j) that depend on are much 
smaller than 1. That is why means (R 2 ),^ and (R 2 ) (j ^, which are expressed via {Nk)/^ by 



eq. (|56j) . depend only slightly on z^. In this case the variances 

Var (x + y) « 2 7 " x |zi| 2 , Va ri (x + y) w 27- 1 |z 2 | 2 



are much bigger than squares of the corresponding means 



- iey) 



U) 



Let us consider uncertainty relations in the semiclassical ICS. Let F\ and P 2 be two 



self-adjoint operators satisfying the commutation relation 



Pi, P 2 



1F3, where P3 is a 



symmetric operator with a real mean (P3). In this case the uncertainty relation 

1- 



Var(F 1 )Var(F 2 ) > -(P 3 ) 
holds, see e.g. [30]. Adopting this general relation to our particular cases, we obtain: 



Varj (Pi) Vaxj (x + y)> fr 2 (Pi + ieP 2 ) (j) , Vaij (L z ) Vaxj (x + y)>— (x - iey) 



(i) 



(79) 



Here (x — is given by d74l) and (Pi + ieP 2 ) 



(J) 



a.s 



can be represented by the help of (125)) 
2 



(Pi+ieP 2 ) (i) =2 7 ^ (ai) (j) 
The variances Var,,- (PjJ an d Var,, (L z ) can be expressed via the Var^ (Nk) (|7T|) . 

Var, (Pi) = (2 7 /i 2 ) 2 Var, (iVi) , Var, (L z ) = tfVarj (JV X - iV a ) , 
and Var., (a; + y) are given by (|77)l . 



We note that 



that 



z\ - z 2 \ 



(«i) 



a) 



|zi| for any ||zi| — |z2|| , and for defmiteness sake, we suppose 
|z 2 || for ||zi| — |z 2 || <C 1. Then, using (fTTj) and (f77| . we see that for 



I \Zi I — I £2 1 1 ^ 1 the products of the variances from (|79|) are close to their possible minimal 
values, 



Vaxj (Pi) Var 3 (x + y) « 4ft 2 (Pi + «eP 2 ) (i) , Vaxj (L z ) Var 3 (x + y) « ft 2 (x - iey) 



and for ||zi| — \z%\\ <C 1 these variances are much bigger than the means H (Pi + ieP 2 " 
2 



0') 



and (ft 2 /4) (x — ity)ij) , respectively. 

Of course, the AB effect is global. However, there exists a difference how this effect man- 
ifests itself in the pure AB field and in the MSF. In the latter case, there exists a possibility 
to characterize especially constructed quantum states with respect to their "closeness" to 
the AB solenoid. Namely the CS have, in a sense, such characteristics. The more close are 
such states to the solenoid, the more they are affected by it. 
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5 Time dependent CS 



On the base of ICS discussed above, one can construct already time-dependent CS (we 
call them simply CS in what follows) as solutions of the corresponding nonstationary wave 
equations. One ought to mention that CS for nonrelativistic spinless particle in the MSF were 
constructed in our recent work [27] . Below, we are going to construct CS for nonrelativistic 
and relativistic spinning particles in 2 + 1 and 3 + 1 dim. 



5.1 Nonrelativistic particles 

In 2 + 1 dim. the quantum behavior of nonrelativistic spin up particle ( ant ip article) is 
governed by the Pauli equation (|1 1 1(1 . where the Hamiltonian can be represented as follows 



N t + (1 - ae) /2 



<T = ±1 



Solutions ^5a(t, r) of such an equation are 

*2 P (i, r)=Afex.p{-i[oj NK (a-e) /2]i}$ ff (t, ip, p)v a , 

respectively for ± cases, where Af is normalization constant, v a is given by (|99l) . and func- 
tions $ CT are solutions of the following equation: 

idt<S>a{t, (fi, p) = <TOJNB.Nx<f><r(t, If, p) . (80) 

One can obey (150"|) setting <& CT (t, <p, p) = &zi, z 2 ,a( l P, p) , where Z\ (t) is a complex 

Zl=Zl(t) 

function of time t. Then 

id t^zuz 2 ,A^ P) = i ^ 9 ^z l ]z 2 ,a(^ P)> h = dz 1 /dt. (81) 



Substituting (|8T|) into (|80j) . we find iz\ = ctwnr^i, where (|43[) is used. It is convenient to 
write a solution for Zi(t) as follows: 

zi(t) = -\zi\exp(-i<np), ip = w NR t + tp , (82) 
where \zi\ is a given constant. Thus the functions 

*& P (*. r) =^exp{-i[u; NR ( ( r-e)/2]i}<f (t)!Z2)ff (^pK (83) 

are solutions of the 2 + 1 Pauli equation for spin up particle. At the same time they have 
special properties that allow us to treat them as CS and even SS under certain conditions. 

Consider the 2 + 1 Pauli equation (jlllj) for spin down particles. The corresponding 
Hamiltonian reads: 

Ni + (1- at) /2 



Solutions vj/d° wn (i ; r) of such an equation have the form: 



#down (t) r) = ^ ex p{-i[a; N R (a + e) /2]t) $_ CT (<, <p, p)v_ a , 

where a = +1 for particle and a — — 1 for antiparticle. Similar to spin up case, one can 
construct CS as follows 

* ( &T n & r ) =^oxp{^[u NR (a + e)/2]t}^ i \ t) z2 ^, p)v_ a . (84) 
We note that means (F)^ in such CS are reduced to given by (|47p. 
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In (3 + l)-dim., one can find CS for nonrelativistic spinning particles, with a given 
spin polarization s. Such CS obey the nonrelativistic Dirac equation with the Hamiltonian, 

(a±P±) 2 + P| /2M, see the Appendix A, and have the form: 



} *c& r), 

where representation (157)) is used. If to consider only physical observables F that do not 
depend on z, then means (-P 1 )^-) in CS (|85[) are expressed via the corresponding means on 
the xy-plane, i.e., via the corresponding means (|47p for (2 + l)-dim. particles. 



(j')NR 
CS, er.s 



(x) = exp 



(gg) crt 
2M 



5.2 Relativistic particles in (3 + l)-dimensions 

In 3 + 1 dim., we consider the Dirac equation in the light-cone variables, such that x- plays 
the role of time. Solutions VP of such an equation have the form: 



*\A X ) = A/"exp<j - — 



/Mc 



fiw (1 — ae) I a;_ 



-av a 



where 



x $\ t<T (x- t (p, p) 

.d$\,<j(x-,<p, p) 



dx~ 



G)Ni$\, a {x-,<p, p) , 



see the Appendix A. 

One can obey ((86)) setting 



Zx=Zl(X-) 



(86) 

(87) 
taking into 



where z\ (xS) is a complex function of the time X-. Substituting (J57J into 
account 

gggU^y, g) y) 

and (I43p . we find idz\/dx- = Cjz\. A solution of the latter equation has the form: 

zi (x_) = — \z\\ cxp(— iip), ip = Qx- + i/jq, (88) 

where \z±\ and ijjo are assumed to be some constants. Thus, we have a set of solutions of the 
Dirac equation in the following form 



*cs,<» =A/"exp 



2ft 



Mc 



+ ftw (1 — ere) I 



x $ 



0) 

Zl(x_), Z 2 



(89) 



We interpret these solutions as CS with light-cone time X- evolution. 

Suppose, we deal with physical observables F that do not depend of x+, which is nat- 
ural for the axial symmetry of the problem under consideration. Matrix elements of such 
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observables in CS (|89)) (we use the inner product (|126[) on the hypersurface x_ = const) 
take the form: 

where the inner product (, ) ± is given by eq. (|105[) . That is why means (F)^ hi such CS 

are expressed via (F)(j\ given by (jT71) . 

Following the same way in the spinless case, one can construct CS that are solutions of 
the Klein-Gordon equation. 



5.3 t and rc_ evolution of mean values 

Let us calculate means {x)^ and (y)rj\ hi the nonrelativistic CS constructed above. These 

means are expressed via the means (x — iey),^, which have the form (|53|) . (|54|) . Taking into 

account eq. (1821 . one can see that means (x)^ and are moving along a circle on the 

xy-plane with the cyclotron frequency wnr, he., the trajectory of the means has the classical 
form. The same equations allows one to find a mean radius (R)^ of such a circle and the 

distance {Rc)(j\ between its center and the origin, 



(i?) (0) = v/2^NAi-^(M,M), (Rc\ 0) = var^N. 

CR) (1) = v^M, = V2r :T N A Ma (|z 2 |,N). (90) 

Note that for the spinless particle \i a = /i. 

2 2 

In the general case, the quantities (R)(j) aid (R c )^ do not coincide with the correspond- 
ing quantities (R 2 )(j) an d (R%)tj) given by eq. (|56|) . The latter quantities are expressed in 
terms of means of square of the transverse kinetic energy and J z according to (l2l?|) and P5)) . 

It follows from eq. (|54|) that Ai_ MtT (|zi|, |z 2 |) < 1 and A^d^l, |zi|) < 1. This allows us 
to give the following interpretation for two types of states with j = 0,1. States with j = 1 
correspond to orbits that embrace the AB solenoid (which corresponds to \z±\ > hi the 
semiclassical limit). For such orbits (R c )m < Rc, where the quantity R c — \/2h/Muj \z^\ is 
interpreted as a distance between AB solenoid and the orbit center (see the classical limit 
of eq. (H2)). At the same time, the mean radius of the orbits coincides with the classical 
radius R — W2h/Mu> |zi|. The interpretation of R as the classical radius follows from eq. 
((62)) in the classical limit. States with j = correspond to orbits that do not embraces the 
AB solenoid (which corresponds to \z\\ < \z%\ in the semiclassical limit). For such orbits 

(Rc)(o) = Rc and ( R )(o) < R - 

One can see the standard deviations Sj (R), 5j (i? c ), and Sj (x + y) = iJYaXj (x + y) in 
CS (f72l and (177)) . are relatively small for the semiclassical orbits situated far enough from the 
solenoid, i.e., for ||zi| — \z%\\ 3> 1. In this case the CS are in main concentrated near classical 
orbits. In the most interesting case when a semiclassical orbit is situated near the solenoid, 
such that the condition ||zi| — \z2\\ -C 1 holds, the standard deviation Sj (x + y) increases 
significantly, Sj (x + y) = 5' (R) f=s 2n~ 1 / 4 j~ 1 / 2 \zi I 1 / 2 , while the standard deviations Sj (R) 
and Sj (R c ) remain relatively small. In this case R ~ R c , however, (R c )^ < R and (R)r \ < 
R c , as it has to be for the semiclassical orbits. Thus, the standard deviation Sj (x + y) of 
particle positions near classical orbits is relatively large at R ps R c , such that 8' (R) 3> 

R — (R c ) (1) , (R)(q) — R c ■ We show the corresponding spreads on Fig. [2] (where R c = 



(R c ) (1) andi?=0R) (o) ) 

We stress that for fi ^ 0, relations between CS/SS parameters of particle trajectory in 
a constant uniform magnetic field differ from classical ones due to the presence of the AB 



22 



Figure 2: Spreads of particle positions near classical orbits for R ss R c . 

solenoid. Above, we have demonstrated this, considering the radius R (related to the energy 
of particle rotation) and the distance R c (related to particle angular momentum). Such 
relations do not feel the presence of AB solenoid for fx — 0, and, even for /i ^ 0, in the 
classical limit. 

For relativistic particles in (3 + l)-dim., we consider means (x — iey)^ in CS (|89p on 
the hypersurface X- = const. Such means are reduced to the means (x — *ey)(j) represented 
above by expressions (l53l) . (l54l) . The relations (l90l) and ((56)) remain true. Here however, the 
evolution is parametrized by the light-cone time x_, via the function zi(x-) given by eq. 
(|88j) . One can see that means {x)^ and (y)^ are rotating along circles on the xy-plane 
with the synchrotron frequency oj, i.e., their trajectories have the classical form (|14[) . 

5.4 Ultrarelativistic particles in (2 + l)-dimensions 

In subsec. 15.21 we have constructed relativistic CS in 3+1 dim., and in subsec. 15. 3[ we have 
demonstrated that in such CS the means have the classical form (fTi)) . We succeeded to do 
this using light-cone parametrization of the evolution via the function zi(x-) given by eq. 
(|88j) . Such a parametrization is possible only in the relativistic case in 3 + 1 dim. Indeed, 
using (|131|) . we can represent eq. (|122[) for eigenfunctions of Po + P3 with the eigenvalues 
A in the form of the first order Schrodinger-like equation (|132j) . where x- plays the role 
of the time and the operator H x _ — Q 2 (2\Mc)~ plays the role of the Hamiltonian. The 
Hamiltonian H x is quadratic with respect to the momentum operators. In the cases in 
(2 + l)-dim. considered above, the light-cone variables x± cannot be introduced. Then we 
have to use the time t parametrization of the evolution. This is the reason why we cannot 
construct CS as exact solutions of the Dirac equation. This is a consequence of the fact that 
in the case under consideration, Dirac Hamiltonian is not quadratic in momenta and the 
corresponding ICS do not maintain their form in course of the evolution. Below, we consider 
an example of such an evolution of ICS. We take massless £ = +1 fermions in (2 + l)-dim. 
with the Hamiltonian H § — ccrP±. 

One can see (using results of the Appendix A) that ICS ([38]) obey the following relation 

fa P) = ±cfi *g+ 1 , ) ft (V, P), (91) 
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where IIo = IIo (^1) i s given by ([3"5"f at a — Then a formal solution of the Dirac 

equation, with ICS as an initial condition, reads: 



±,z 



i.zafoP) P) = ex P 



±, Z1,Z2 



(v>, P)- 



(92) 



We call such solution quasi-CS in what follows. As usual, we define means of an operator 
F in quasi-CS flU]) by W)(j)> 



(F(t)) 



OKI w 



a,±) 



\ ±, Zl ,Z2 * ±, ^1 ,#2 

One can see that in the semiclassical limit the means 



(x (t) - iey (t)) . ±) = (oa (t))y, ±) - (oi (t)) & . ±) 



(93) 



(94) 



are moving along the corresponding classical trajectories. Since the operator a-i commutes 
with Ni, the mean (ci2 -m does not depend on time and coincides with its initial value, 

( a 2)(j) ' the latter is given by eq. (1551) . Calculating the mean (a\ (t)) (A _1_ \ j WG 



(°2 W)(i,±) 

find that 



(ai (*)) 



the latter is given by eq. (|55|) . Calculating the mean (a\ (i))(j ±y 



(0,±) 



(«1 (*))(! 



Ea=±lQl-^(l Z lU Z 2|) 

£ CT=±1 exp [±»n (^^) t] Q^(N, 



(i,±) 



E CT =±lO/i„(N,N) 



(95) 



see the Appendix B, where the frequency operator fl {Nij is given by eq. (|146l) . 

Let now \zi\ 2 ^> 1, which corresponds to the semiclassical limit. Using relations (|59"|) and 



60]) . we represent (|95[) in the form 



(«i (*)) 



(o,±) 



(Ol (*))(i ±) = 2 



= M«i 

The semiclassical expansions of the operator (7 is: 

1 d 



E„= 


,±i exp 






1 


Ea=±lQl-^(l 2; lU 2; 2|) 


E CT= 

* 


,±i exp 






i 


E CT =±iO/^(M,M) 


N 


d V 






2 d 









f2 ss w (Ull 



1 



4|«i|d|«i| 



ON 



, W (|^|) = c|gB|f- 1 (|z 1 |), f (1^1) = V2fic| 9 B| 



such that the standard deviation of f2 is of the order w (|zi|) \ z%\ . Here £ (|zi|) is the mean 
energy in quasi-CS, in the classical limit. Taking into account next-to-leading corrections to 

and using decompositions (jMJ) and (1551) . we obtain (ai ±) = ( a i)(j) e±? '^ i * i 



exp 



where 



Oj =w(\zi\) [l + O (\ Zl \ 2 )] , ||zi| - \z 2 \\ > 1; 
(-l) j 



cjj = a; (|^i |) 



o(W\- 



21 - M < 1. 



(96) 
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In such an approximation, lj (\zi\) coincides with the classic synchrotron frequency u given 

by eq. ©. 

Thus, in the classical limit, means {x)(j _y and {y)(j ±\ in the quasi-CS are rotating along 

circles on the xy-plane with the synchrotron frequency w. The greater is |zi| 2 the smaller is 
spreading of the mean trajectories. One can see that additional (due to the next corrections) 
modifications of quasi-CS are essential only for evolution time that is much greater than the 
classical rotation period. 

The mean radius of a trajectory is 



(ip; 



2(^)^+1+06 



7 -i. 



The next-to- leading corrections to this relation can be find by the help of eqs. (|51|). (p5)) . 
and (|67|) . We note that in each approximation the classical relation between the rotation 



frequency and the radius holds, such that ujjJ {R 2 )(j ±) = c (we remind that in the graphene 
case c means the effective velocity of light, that is the Fermi velocity vp). 



6 Summary and Discussion 

A new approach to constructing CS/SS in MSF is proposed. The main idea is based on 
the fact that the AB solenoid breaks the translational symmetry in the xy-plane, this has 
a topological effect such that there appear two types of circular trajectories which embrace 
and do not embrace the solenoid. Due to this fact, one has to construct two different kinds 
of CS/SS, which correspond to such trajectories in the semiclassical limit. Following this 
idea, we construct CS in two steps, first the instantaneous CS (ICS) and the time dependent 
CS/SS as an evolution of the ICS. 

The approach is realized for nonrelativistic and relativistic spinning particles that allows 
us to build CS both in (2 + 1)- and (3 + l)-dim., using some universal constructions, and 
gives a non-trivial example of SS/CS for systems with a nonquadratic Hamiltonian. 

It is stressed that CS depending on their parameters (quantum numbers) describe both 
pure quantum and semiclassical states. An analysis is represented that classifies parameters 
of the CS in such respect. Such a classification is used for the semiclassical decompositions 
of various physical quantities. 

In the pure quantum case, the mean values depend significantly on the particle spin 
and on the mantissa and are quite different from the corresponding classical values. In 
the semiclassical approximation, relations between CS/SS parameters and parameters that 
characterize classical trajectories are established. In the general case these relations differ 
from ones in the pure magnetic field and such a distinction can be treated as AB effect in 
the CS/SS. The classical relations correspond to the leading approximation for sufficiently 
large radii. Thus, the leading approximation in the semiclassical expansions corresponds 
to the classical limit. Next-to-leading terms define physical quantities in the semiclassical 
approximation. These terms depend on the space dimension and particle spin. 

The following properties of the constructed time dependent CS/SS should be stressed: 

a) In the nonrelativistic case, both in (2 + 1)- and (3 + l)-dim., the time dependent 
CS in each time instant retain the form of the corresponding ICS. The mean trajectories 
in such CS coincide with classical ones, whereas the particle distributions are concentrated 
near the classical trajectories in the semiclassical approximation. In the presence of the AB 
solenoid, the spread of particle positions near the classical trajectory depends essentially on 
the mutual disposition between the trajectory and the solenoid. Such a spread is growing 
for trajectories situated near the AB solenoid. It should be note that namely due to the 
bounded character of particle motion in the MSF, particle positions are essentially sensitive 
to the topological effect of breaking of the translational symmetry in the ir?/-planc due to 
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the presence of AB solenoid. Thus, in spite of the well-known fact that the AB effect is 
global, in the MSF quantum states can be classified with respect to their " closeness" to the 
AB solenoid. 

b) In the relativistic case, in (3 + l)-dim., CS are constructed in the light-cone variables, 
where the evolution is parametrized by the light-cone time X—. Such time-dependent CS 
obey all the properties as CS from the previous item a). 

c) In (2 + l)-dim., we constructed time-dependent SS for massless fermions. Such a 
problem can be related to the graphene physics. We call the constructed SS quasi-CS since 
they retain ICS form with time evolution in the next-to-leading semiclassical approximation. 
In such an approximation, the classical relation between the rotation frequency and the 
radius holds, the rotation frequency coincides with the classic synchrotron frequency in the 
leading approximation. We stress a principle difference between (2+l)-dim. and (3+l)-dim., 
in (2 + l)-dim. Dirac Hamiltonian is not quadratic in the momenta, whereas in (3 + l)-dim. 
it is. Namely this fact is responsible for the destruction of ICS in course of the evolution. 
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Appendix A. Quantum stationary states 

.0.1 (2 + l)-dimensions 

In (2 + l)-dim., the total angular momentum operator J = —ihd^ + ha 3 /2, which is a 
dimensional reduction of the operator J z in (3 + l)-dim. (z-component of the total angular 
momentum operator given by (|22[)). is self- adjoint on Dfj and commutes with . There 

exist common eigenvectors ipni, n 2 {f: p) of operators H ® and </, 

H^ln 2 fa P) = c Po ^l n2 (<p, p), c Po =±£,£= \J (Mc 2 ) 2 + £ 2 , 
■W, n 2 fa P) = J^l n 2 fa p),J = eh{l -l + 1/2) , j = 1, 2. (97) 
It is convenient to use the following representation 

^InM P) = [°* fa - ^i) + Mc] U^fap), 

^InMP) = E C ^ { nln 2 ,A^ PK, (98) 
a=±l 

where 

«*-(£). *-!=(?)• <"> 

(i) 

and Co- are some constants. The columns Wn lj n 2 are solutions of the eigenvalue problem 

c 2 (aP ± ) 2 n» n2 fa P ) = Slu^l „ 2 fa p). (100) 

We note that the relation (crP^J = — ehc^ 1 \qB\ <r 3 holds, which gives a relation to 
the energy spectrum for the spinless case. 



2G 



The functions &n}, n 2 ,a-( ( P, p) have the form 

P) = -^exp{ie p - I + (1 - ea) /2] <p}/ n2) m (p) , 
ni=m, n 2 = m- 1 + (1 - ecr) /2 - p, 1 = I - (1 + e) (1 + cr) /2, -oo < Z~ < - (1 + tfe) /2; 
^U,^' P) = AAexp{ie [/ - / + (1 - e<r) /2] <p + ien [I - (1 - e) (1 + cr) /4]} J m , n2 (p) , 
ni=m + f-(l-eo-)/2 + /i, n 2 = m, (1 - i?e) /2 < I < oo ; p = |r 2 . (101) 

Here I n ,m(p), m > are Laguerre functions that are related to the Laguerre polynomials 
L^(p) (see [21]) as follows: 



i m+a , m ( P ) = / r ( m + t+i) e " P/V/2£ - (p) ' ^ = ^1^^ ' (102) 

and A/" are normalization constants. For any real a > — 1, the functions I a +m,m(p) form a 
complete orthonormal set on the semiaxis p 0, 

/*°° °° 

/ Ia+k,k{p)I<x+m,m{ft)dp = h,m , /] Ig+m,m(p)Ig+m,m(p') = 5(p ~ p') ■ (103) 
Jo m=0 

Thus, the domains Dfj are described completely by asymptotic behavior of the functions 
from eqs. (] 1 1 1) . 

We define the functions $^f?+ ai „ 2 + S2 CT associated with transformation (|26|) as follows: 



(0) 

rii+si , ri 2 +s 2 



(1) 

ni+si , n 2 +s 2 



a ((p, p) = Afexp{ie [l - I - Si + s 2 + (1 - ea) /2] <p} 7„ 2+S2i , il+Sl (p) . 
CT (cp, p) = A/"exp{ie [Z - I - s x + s 2 + (1 - ecr) /2] ip 



+7T [I + ai - s 2 - (1 - e) (1 + a) /4]} 7„ 1+Sl ,„ 2+S2 (p) , s x = 0, ±1, s 2 = 0, ±1. (104) 

There appear new functions $^ _ x CT with n 2 = m + 1 + (i9 — cr) e/2 — p and [ _ 1 n2 g 
with ni = r7i + (a — d) e/2 + p. Such functions were not defined by eqs. (jlOip . In addition, 
for ri\ = or 7i 2 = 0, one has to bear in mind that 



ai 



m=0 



0. 



(0) 



This allows as to interpret $>n{ 



and $ 



(i) 



m=0 

as vacuum states. 



Let us define an inner product of two functions /(<p, p) andg((p, p) as 

^ ,9 ^- L = 2nJ [ dp J dlfi f*&'P)9&'P)- 
With respect to such an inner product, the functions (|101[) form an orthogonal set, 



(105) 



(106) 



Moreover, these functions form a complete orthogonal set in L 2 (R 2 ). 

The energy spectrum of the self adjoint Hamiltonians (f2"Tj) can be represented with the 
help of the eigenvalues £\ as 



£ 2 ±M =2hc\qB\[ ni + (l-*e)/2}, 



(107) 



where, depending of the t? and e, the quantum number ni takes its possible values according 
to (|101l) . In the general case, eigenvalues £ of the Hamiltonian are expressed via £]_ for 
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cr = +lor<7 = — 1, according to eq. (|97|) . that is why are labeled by the subscript a. 

Irregular at the origin radial functions appear in the domain for I = and a = — 1, and 
in the domain DJj for / = and a = +1. 

We note that depending on ?9e, energy levels of states having irregular radial functions 
coincide or differ from the Landau levels. In any case the difference always depends on fi 
only. Energies of states with j = 1 differ from the Landau levels, whereas energies with j = 
coincide with the Landau levels. If £±/ a ) ^ 0' the com plete set of eigenvectors ipni, n 2 is given 
by eq. (j9"8")l . where constants c CT are arbitrary, e.g., either c+i 7^ or c_i 7^ 0. However, if 
£j_r a ) — 0, a completeness of the eigenvectors implies a special choice of c CT , namely: c_i 7^ 
if e = — 1 and c+i 7^ if e = +1. In this case the only negative energy solutions (antiparticlcs 
with cpo = —M) are possible if e = — 1 and only positive energy solutions (particles with 
cpo = +M) are possible if e = +1. They coincide with the corresponding spinors Un{, n 2 {<P, p) 
up to a normalization constant. This is a manifestation of the well-known asymmetry of the 
energy spectrum of 2 + 1 Dirac particles in the uniform magnetic field. We see that the same 
asymmetry takes place in the presence of the AB-field. For particles, we select c + i 7^ 0, 

c_i = 0, then their energy spectrum is £ — >J [Mc 2 ) 2 + and for antiparticles c+i = 0, 



c_! 7^ 0, then their energy spectrum is £ — yj (Mc 2 Y 

If we define the inner product of spinors ijj(ip, p) and ip'((p, p) as follows: 



£ 2 



1 
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then the inner product of eigenvectors ipfii 



dp 



U) 



2tv 



dipip^(<p,p) ip'(if,p), 
has the form 



\i n' x , n' 2 ' "ni 



D 



2Mc U 2 (ap + Mc) (sj/^, ^?, n2)(T 



(108) 



(109) 



where the inner product (, ) ± is given by eq. (|106[) . With respect to the introduced inner 
product, eigenvectors (|98|) form an orthogonal set for any 

By the help of eigenvectors (|98|). we obtain the following solutions of the Dirac equation 
with a given energy cpo = ±£, in (2 + l)-dim.: 



po . m, I 



exp 



■ (cpot) 



^n>. P)- 



(110) 



We believe that nonrelativistic motion is described by the corresponding Pauli equation 
in (2 + l)-dim., 



/2M. 



(Ill) 



Solutions of such an equation can be obtained from (|110[) in the nonrelativistic limit. They 
have only one component, taking into account that in (2 + l)-dim. there is only one spin 
polarization. Let us consider, for example, spin up particles (C = +1). Then we obtain from 
(EH): 



*? >P Ki,r) 



exp 



4> 



(i) 

n 1 , n 2 ,± 



i(¥>, p)v±i, 



(112) 



where £, 



(±) = ^J-(±i)/^^ c2 - The corresponding inner product of these solutions reads 



(j')"P 
±, m', /' 



(j)up _ 



D 



(./) 



ni, rijjil ' Til, ri2,±l 



(113) 



By the help of relation (|20j) . we obtain solutions that describe spin down particles: 



^(J)down 



exp 



NR j 



0) 



^(j')down ^(j)dowiA = (&W) $ 0") A 

V ±,m',i' ' ±,m,( / n \ »»i)«2.Tl' nx,n-i,^\ ) 



(114) 
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It is worthwhile to make the following remark: Usually, in the nonrelativistic limit, 
negative energy solutions (those which after the charge conjugation operation represent 
wave functions of particles with opposite charge —q, and which are, in fact, antiparticles) 
of the Dirac equation are not considered. It is supposed that all the information about the 
quantum motion of the antiparticles can be extracted from particle motion. The latter is not 
true in the case under consideration (for Dirac eq. with MSF in (2 + l)-dim.). Here energy 
spectra for particles and antiparticles are quite different. This asymmetry was the reason 
for us to represent explicitly, even in the nonrelativistic limit, the negative energy solutions 



(j)up 



At, r) and * 



(j)down 



(i,r), which correspond to spin up and down antiparticles. 



- , m, I —,m,l 

For massless £ = +1 particles and antiparticles (in what follows we call all such particles 
simply fermions) in (2 + l)-dim., the self-adjoint Dirac Hamiltonian is H' = ccrP±. Its 
eigenvalues are cp = ±£ , where £ = (the eigenvalues £u<r) are given by eq. (|107[) ). 

The corresponding eigenvectors of H® have the form: 



^+%i,r) = exp 



T- h {Et) 



.O'.+i) 

±,ni, n: 



(<p>p)> 



,(0, + l) 



m'miViP) 



(0) 



±^i 1 ) + (l- e )/2,n 2 ,-l(^' P) 
(1) 



u ±,ni,n a (V.p) = . _(i) , J. £ r o. 



where &n},n 2 ,o- are given by eqs. (1101 1) . The inner product of such solutions reads 



\ ±,m',i" ±,m,lj D V »i+(H')(l 



0) 



+cre)/2,ra^,cr' ni+(l— j)(l+cre)/2,n.2,CT 



tr=±l 



(115) 



(116) 



In addition, there are nontrivial zero-mode (£ = 0) solutions: 



(0,+l) (0,+l) 
U 0,n' n' ' M 0,ni,n 2 



V «1>™2 



We , 

m=0 
$(°) 



m=0 



(117) 



As follows from (|2"D|) , for massless £ = — 1 particles the corresponding eigenvectors can 
be represented as 



Vjt.r) = exp 

,(o,-i) 



(118) 



The inner products of these solutions coincide with ones of solutions given by eq. 

■mm-. 



\ ±,m',V ±,m,lJ D \ ±,mV ±,m,lj D ' \ 0,n' 1} 



(V°< +1 



) ..(0. + 1 ) 



D 



.0.2 (3 + l)-dimensions 

Here we consider Dirac equation (|16[) in (3 + l)-dim. Let us introduce projection operators 
P(±) and two kinds of Dirac bispinor ^(±), 

V {±) = (1 ± a 3 ) /2, (P (±) )t = p (±) , (7> (±) ) 2 = V {±) , V {+) V ( -) = 0, V {+) + V { -) = I, 
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where I is unit 4x4 matrix, such that any \1/ can be represented as = ^(+) + V E'(_), ^(±) — 
P(±) 1 ^- Then Dirac equation (fTBT) is reduced to the following set of equations 

(p + = (A - - Q*(+); 

Q = (a x Pi) + Mc7°, P = iKDo, (119) 

where a 1 = 7 7*. Due to the axial symmetry of the problem, it is convenient to use the 
following representation for 7-matrices (see [9]), 

7 = diag (a 3 , -a 3 ) , 7 1 = diag (ia 2 , -ia 2 ) , j 2 = diag (—icr 1 , ia 1 ) , 7 s = antidiag (— J, 1) , 

(120) 

where / is unit 2x2 matrix. Nevertheless, expressions for a 3 and 7 5 are the same in 
the representation (|120[) and in the standard representation. 
One can see that 

Q 2 = M 2 c 2 + Ql, Q 2 ± = (a ± P ± j 2 = Pi - ehc- 1 \qB\ E z , (121) 

where Y, z — diag (a 3 , a 3 ). In the MSF, the operators Po + P3, Po — P3, and Q mutually 
commute, such that (as it follows from (|119l0 bispinors &r±) obey the same equations: 

(Po 2 -P 3 2 -Q 2 )*( ± )(z)=0. (122) 
Representing ^(±) via spinors ip and x, 

*(-) - ^ ( -% ) . = K A ) ' (123) 

we find the following equations for the spinors: 



P 2 -P 2 -(Pl-e\ qB \^ 3 )~M 2 c 2 



= 0, (124) 



(P + P 3 ) X (x) = a 3 <rP ± + Mc) ij,(x). (125) 



We note that both solutions and V E'(+) enter into a complete set of functions on the 

hypersurface t — const. 

The inner product of Dirac bispinors on the light-cone hypersurface X- = const has the 
form: 

(*,#')*_ - J ^V ( ^'dx+dx l dx 2 = ^ J (*(_), *'(_)) (126) 

see |29j , where the inner product of four-component spinors ^ and *f/' on xy-plane is defined 
as 

= — J dpj^ tkp&fapl&fap). (127) 

It is expressed only in terms of the components At the same time, a complete set of 

functions on the hypersurface X- — const consists only of ^(-V 

In the case under consideration, the operators Po, P3, J z = L z + S z /2, and a spin 
operator S z (z-component of a polarization pseudovector) 

S z = i (P"E Z + E Z P") = 7°S 2 Mc 2 - 7 5 cP 3 , (128) 

are mutually commuting integrals of motion (all these operators commute with the Hamil- 
tonian H ) [Tj2 [16]. In addition, the set Po, P3, J 2 , £ 2 , and Qj_ represents mutually 
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commuting operators, which, at the same time, commute with a 3 . This fact allows one to 
find solutions that are eigenvectors for the latter set. To this end one has to subject 
spinors ip to the following equations: 



(Po + P 3 ) i>{x) = \MciP(x), J z iP(x) = J z i/j(x), 
c 2 (aP^%{x) =£l [rT) ^{x), a^{x) = aiP(x), 
where £±r a ) is given by (|107l) . Thus, we obtain for ^(_) : 

(Po + P 3 ) * ( _) = AAfc* ( _ ); J,* ( _) = J^ ( _), J z = eh{h - I + 1/2), 



(129) 



£V)*<-), 



(-)' 



±1. 



In the light cone variables (fT3")) . we have the following representation 



P + P 3 = 2ih 



d 

dx+ ' 



Pn - P, = 2ifr 



_9_ 

dx - 



(130) 



(131) 



Then, we can represent eq. (|122[) for eigenfunctions of Po + P 3 with the eigenvalues A in the 
form of the first order Schrodingcr-like equation 



2ih\Mc- 



d 
dx- 



*(_)A(a:)=Q, 



(132) 



and we find a complete set of solutions \P V_) in the following form: 



(./) 



( — ) A, m, Z,er 



\Mcx+ 



(Mc 



hu (1 — ae] 



X<P, p) 



— aVa 



(133) 



where ui is given by (JTJJ), w CT by eqs. (JMJ, $n?,n 2 ,<T by (|101[) . and A > for particles 
and A < for antiparticles. We note that the quantum number A is associated with the 
corresponding classical quantity A from (|12|) . 

We note that the spin integral of motion S z does not commute with a 3 such that solutions 
(|133|) are not eigenvectors of S z . One can use the operator T, z instead of S z to characterize 
the spin polarization. In spite of the fact that S Z ,P 7^ and, therefore, S 2 is not an 
integral of motion with respect to the i-evolution, E z is an integral of motion with respect of 
the evolution in the light-cone "time" X-. That is why the "spin polarization" a of solutions 
(|133[) is conserved with the time X- . Taking all this into account, one can calculate the light- 
cone inner product (|126[) of solutions (|133[) : 



V (-)A',m',/',(T' : 



(-)A, 



(An) 2 h 



^(A'-A)(^P n , iff ,$« n2 , ff )^, (134) 



where the inner product (, )j_ is given by eq. (|106[) . 

Let us consider the quantum motion of spinning particles in the nonrelativistic limit. 
To this end it is more convenient, instead of solutions (|133l) . to use another set of solution 




-(cp t + p 3 z) 
n 



* s (x x ), s = ±l, 



1+ [p 3 /c + sM)/M ip Po<s (x±) 



-1 + p7c + sM /M tp P0 , s (x ± ) 



(135) 
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where M = \/ M 2 + {ps) 2 ■ These solutions are eigenvectors of mutually commuting integrals 
of motion Pq, P3, J z , and S z , 



P ^s(x) =po^ s (x), P 3 V s (x) = p 3 ^ s (x), 
JA(i) = ^f,W, Jz = eh(lo — I + 1/2), S 2 f s = S Mc 2 * s , 
The spinors ippo,s{x±) obey the equation 

<tP_l + sMcer 3 ) Vwfa-O = Po^ Po ,s(x±). 



(136) 



(137) 

One can see that at fixed s and p 3 , eq. (I13T[) is similar to eq. (TT51) in (2 + l)-dim. such that 
its solutions will be used in what follows. 

In the nonrelativistic limit, the spin operator S z is reduced to Sf R = 7°E z Mc 2 and 
71/ = M. Then for s = +1 eq. (|137p coincides with eq. (TT51) . We remark that ifj p0t -i(x±) — 
a i/)- POt i(x±) . As a result, we obtain wave functions of nonrelativistic spinning particles 
from eq. (|135[) : 



C?')NR 
±,]?3,m, l,s 



, (x) = ex.pl — - 



, (Pa) 2 * , 
± — +P3Z 



2M 



0")NR 
±,m, l,+ 



7 (*, r) 



* u;uiJ r(t,r) 



0")NR. 





w Cj)down ( t v 



(138) 



where spinors 

*2 } ^ P r and * 

are respectively solution (II 12|) and (|114[) of the Pauli 

equation in (2 + l)-dimensions with MSF. Thus, wave functions of nonrelativistic spinning 
particles (antiparticles) in (3 + l)-dimensions obey the nonrelativistic Dirac equation with 

the Hamiltonian = (q\ + P|) /2M. 



Appendix B. Mean (a\ (t))^ ± \ 

To study the mean (ai ±) > one nas ^° calculate the matrix element y&±' t ^ Z3 (t) , o,\^±'^l 

The latter can be reduced to a matrix element with respect to the initial ICS 2:1 22 (y> z 9 ) 
as follows: 



where 





ic s 




ic - 


exp 


±-n t 

n 


a\ exp 


n 



The operator a{ (f) obeys the equation 

da\ (t) 



Let us consider the commutator 



dt 



a\(t),tl 



a] (t) , n 

First, we write 



a\ (t) , fto| = [at (t) , ngl n„ 1 + n% \a\ (t) , n^ 1 



(139) 
(140) 

(141) 



(142) 
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where the identity 



tator 



a}(t), fto 1 



a\(t),tl \ = \a\ (t) ^IIIq 1 
as follows: 



is used. Then, we represent the commu- 



at (t) , fto 1 ] = A jf ( a t (t) _ gj (r2) ) e -ft^ 
Si(r 2 ) =e- ft o- 2 al(t)e+ ft o- 2 . 



The operator 6 J (r 2 ) obeys the equation 



t („2\ 



db\ (r 
dr 2 



&I(r 2 ),n 2 



M ^ (r2) 



and coincides with aj (i) at t 2 = 0. Such a solution of eq. (|144|) reads: 

SJ (r 2 ) = a{ (t) exp (-2ft r 2 /c) . 
Substituting (|145[) into (| 143)) and calculating the integral, we obtain: 

at (i), fto 1 ] L ' ' - x 

Using this result in (jl42[) . we find 



(143) 



(144) 



(145) 



fto 1 - (tl 2 Q + 2h\qB\/ t 



at (t) , fio] = al (t) ^n - yng + 2ft / c 

such that eq. (|141[) has the following solution: 



a\(t) = a\e in ( ffl ) t , SI = 2 |gB| ( ^ftg + 2ft | g B| /c + n 



(146) 



where f2 can be interpreted as the frequency operator. With account taken of (|146[) 

in eqs. (|139[) . and using relations (|4*5|) and (|4"9")l . we finally obtain expression (j9"5j) : 
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